Thermoelectric transport near the pair breaking quantum phase transition out of a 

tf-wave superconductor 



o 
o 

(N 
S3 

a 



i 

1/3 

a 

I 

S3 

o 

(N 
> 

o\ 
in 
o 

in 
o 

a 

S3 
O 
O 



X 

S3 



Daniel Podolsky 1 , Ashvin Vishwanath 1 ' 2 , Joel Moore 1,2 , and Subir Sachdev 3 
1 Department of Physics, University of California, Berkeley, CA 94720 
2 Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 
3 Department of Physics, Harvard University, Cambridge, MA 02138 
(Dated: Printed February 2, 2008) 

We study electric, thermal, and thermoelectric conductivities in the vicinity of a z = 2 
superconductor-diffusive metal transition in two dimensions, both in the high and low frequency 
limits. We find violation of the Wiedemann-Franz law and a dc thermoelectric conductivity a that 
does not vanish at low temperatures, in contrast to Fermi liquids. We introduce a Langevin equation 
formalism to study critical dynamics over a broad region surrounding the quantum critical point. 



INTRODUCTION 

Transport of heat and charge provide a useful probe of 
strongly correlated electronic systems. For example, a vi- 
olation of the Wiedemann-Franz law at low temperatures 
signals non-Fermi liquid physics. One possible origin of 
such a violation is proximity to a quantum critical point 
(QCP), where the presence of low energy critical modes 
leads to novel phenomena. There has been much theoreti- 
cal and experimental attention directed towards electrical 
transport at QCP's, in the context of the superconductor- 
insulator transition in amorphous films, and the Quan- 
tum Hall transitions (see [l[ Q for reviews). Recently, 
experimental investigations of thermal and thermoelec- 
tric properties at quantum phase transitions have been 
pursued with very interesting results. These include mea- 
surements of thermal conductivity in highly underdoped 
cuprate systems near critical doping and of thermo- 

electric transport (Nernst and Seebeck coefficients) near 
a heavy fermion QCP However, there has been rela- 
tively little theoretical work on thermal and thermoelec- 
tric transport at QCP's. In this paper we calculate the 
nontrivial electric, thermal, and thermoelectric conduc- 
tivities ataz = 2,<i = 2 quantum critical point for the 
transition from an unconventional (non s-wave) super- 
conductor to a diffusive metal in 2 dimensions. 

In an applied electric field E and temperature gradient 
VT, the electric (a), thermal (k), and thermoelectric (a) 
conductivities are defined by, 
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where j and jq are the electrical and heat currents, re- 
spectively. Thermal conductivity measurements are al- 
ways carried in open circuit (J = 0) boundary condi- 
tions, so that k = k — aT/a. The difference between 
k and k is negligible in metals, but this need not be 
the case in systems with strong particle-hole symmetry 
breaking, as considered here. Our results are shown in 
Table |U We find that, while the dc thermal conduc- 
tivity is metallic, the dc electric conductivity diverges as 



T — > 0. Hence, there is a strong violation of Wiedemann- 
Franz law (WF) in the bosonic sector, with an apparent 
excess of charge carriers. The total conductivities are 
sums of fermionic and bosonic parts, a — Of + <t$ (and 
similarly for a, Thus, the magnitude of the ex- 

pected WF violation is difficult to estimate. Another 
striking effect is that a$ has a very weak T dependence 
at low temperatures (~ lnlnT), unlike otf, which van- 
ishes linearly with T. Hence, experimental measurement 
of the anomalous thermoelectric conductivity would give 
direct information regarding critical transport properties, 
as well as yielding the dimensionless "thermoelectric fig- 
ure of merit" ZT = a measure of the strength of 
particle-hole breaking. 
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TABLE I: Asymptotic low T behavior of transport coef- 
ficients along the line s = s c in the zero temperature 
(fi/T -> oo) and dc (fi/T -> 0) limits. The prefac- 
tor of a(Q — 0) is given in an extreme low T limit b — 
[tan _1 (l/7?) + 4^/(1 + ?7 2 )]/21riln(A/T). Improved values for 
the fi = results, valid for a broader T range appear below 
Eq. JTU). Explicit expressions for a(T = 0) and k/T(Q = 0) 
are given in Eqs. and ((HJ. 



A QCP with either Lorentz or Galilean invariance will 
have infinite thermal conductivity, because a "boosted" 
thermal distribution will never decay; the same logic ap- 
plies to electric conductivity if there is only one sign of 
charge carrier. We find below that a finite r\ at the z = 2 
QCP regularizes the thermal conductivity. It is also cru- 
cial that the z = 2 theory breaks particle-hole symmetry 
so that the thermoelectric coefficient can be nonzero; the 
z = 2 theory is thus the simplest critical theory for which 
all the transport coefficients in (fT]) are finite. 
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MODEL AND CURRENTS 

Our starting point is an electron model with a pair- 
ing interaction favoring unconventional superconductiv- 
ity, and a disorder potential whose main effect is to render 
electrons diffusive. Then, introducing a Hubbard-Strata- 
novich field $ in the Cooper channel and integrating out 
fermions yields a Ginsburg-Landau action for $,01 







dr / d z r 
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where the dissipative term |<9 T |, shorthand for the Mat- 
subara expression \w n \, arises from the decay of Cooper 
pairs into gapless fermions. Note that, for unconventional 
superconductors, disorder is pair-breaking. This insures 
that the critical theory ([2|) is local[6(. Naive power count- 
ing of model {2} yields dynamical critical exponent z = 2. 
Although disorder dominates extremely close to the tran- 
sition, one can choose microscopic parameters such that 
cq. ([2]) describes a large crossover region near the QCP[|j]. 
For instance, for a clean sample with an elastic mean 
free path of 100 nm, Model J2|) is valid provided T is 
larger than a few millikelvin[15]. This paper focuses on 
transport in this region. Equation ^ also applies to the 
onset of antiferromagnetic order in an itinerant electron 
systemQ, with $ an 0(3) order parameter. Many of our 
results apply to this QCP as well (eg. thermal conduc- 
tivity) but the field $ carries spin, and not charge. 

The electric and heat currents j = and Jq — 
are obtained by making the action ([2]) gauge-covariant 
through the substitution V — > T> = V — ieA e — iAT{idt) 
(see Ref. @), 
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FIG. 1: Phase diagram in the vicinity of the QCP as a func- 
tion of disorder, parameterized by the lifetime r of electrons 
in the diffusive metal. At finite T, interactions shift the 
transition line to the right, so that the diffusive metal phase 
lies above the QCP. The crossover into the quantum critical 
regime occurs when R « T. 



the order of the limits — ► and T — ► 0[8|. Sending 
ri — ^ first yields the dc conductivities, which are more 
readily accessible to experiment, but are typically more 
difficult to compute than their T = analogs. 

For a d = 2, z = 2 theory, the quartic interaction V 
is dangerously irrelevant. Thus, at the T — QCP, the 
correct dynamics is obtained from the limit V — * + . We 
first consider this non-interacting case. Note that, even 
without interactions, finite transport is plausible, due to 
dissipation in action Indeed, a is metallic in the 

zero temperature limit [6(, as shown in the middle column 
of Table |TJ On the other hand, to compute dc conduc- 
tivities, we shall consider finite temperatures, for which 
interactions give important logarithmic corrections. 
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We compute conductivities from Kubo formulas involv- 
ing the dynamical correlations of these currents. 



NON-INTERACTING CASE 



The conductivities are obtained from Kubo formulas, 
G y (Q) = ^j^ n 7 (0 + ifl). For V = 0, n 7 is given by 



a single-loop integral, 



ORDER OFT^O AND -> LIMITS 

Conductivities near the 2d QCP depend on ratios of 
small energy scales (and possibly a UV cutoff scale A), 
and on the dimensionless parameters mV and ij, 
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Here, Go = a, G\ = a, and Gi 
of the external field, and vz = 



: k/T, £1 is the frequency 
1 for the d = 2, z = 2 



model in question. At the QCP, s 



G 7 depend on 



= £E/ H§ {-<Gi^n? (4) 

where Q is the Matsubara Green's functions, Q(k,iv n ) = 
(iuj n — rj\u! n \ — €k — s) , where = fc 2 /2m. Rewrit- 
ing this in terms of the spectral function, Q{k,uj n ) = 

Jtr££^M(fc, - 2 ^.-,, , we obtain, ' 
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In (|5|) we have made the substitution ioj n — > (wi + W2)/2. 
This accounts for the fact that the time derivative in 
jq does not commute with the time-ordering symbol, 
and is necessary to obtain the correct thermal conductiv- 
ity [Io| . Performing the u) n sum and using the identity, 
f°(u> + iVL n ) = f°(oj), before analytically continuing to 
real frequencies, iQ n — > f2 + i5, we find 

im n y (n + id) = -jf ^^--J x (6) 

.A(fc,cj 2 - n)^(fc,wa) (/ Q (w 2 - 0) - /°(w 2 )) , 

where = (exp (/3w) — is the Bose occupation 

factor. The integrals are UV finite, and hence the depen- 
dence on the cutoff A drops out, 



G scaling /on scaling 
7 l"J — S 7 



■ v 



Expression (J6j) can readily be evaluated for generic fi/T. 

Zero temperature (yi/T — ► oo) limit 

The results in the zero temperature limit at the QCP 
(s = s c = 0) are summarized in the middle column of 
Table HI While a is metallic, a and k/T are divergent 
(note that a is purely reactive). The electric conductivity 
depends on the dissipation t), 



(2e) 2 /h rr , 
x | = / ' { [2rj + tt(1 + ij 2 )} arctan rj 

— 7T7/ — r? 2 — (1 + ?7 2 ) arctan 2 77} , (7) 



since rj is marginal. 



DC (fi/T -> 0) limit 



On the other hand, the dc transport properties at V = 



are drastically different. Whereas a = ^x-jpy 

_ 2efci3 1 1 T 



In jjjy are divergent, y is metallic, 
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In order to elucidate the dc results, we employ a simple 
Boltzmann equation approach, which is exact for weak 
dissipation ?y — *■ 0. The dissipating term in the action 
([2]) can be interpreted as the self-energy of bosons due to 
their interactions with the fermion bath, 

T,(iw n ) = r)\u} n \. 



This is purely imaginary for real frequencies (\z\ = 
zsign(Rez)). Hence, while the energy of a quasi parti- 
cle is not renormalized, quasi particles are given a finite 
lifetime. 
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This is the boson lifetime, at bubble level, in the fermion 
model. Here we neglect corrections to the transport life- 
time from higher order diagrams. 

The Boltzmann equation for the occupation function 
/(k, r, t) in the scattering lifetime approximation yields 



dtf + v ■ V r / + k • V k / 



-(/-/°). 



(10) 



The electric and energy currents are expressed in terms 
of /, 



7 Q M) = 
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The conductivities are then 
(2e) 2 -Ul f d 2 k 



(2eK/(k,r,i) (11) 
2 (e k - S K/(k,r,t). (12) 
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For r/ — 0, this gives perfect agreement with the pre- 
vious results for G 7 in the dc limit. Here, we see that 
the reason for the divergences in a and a is that low en- 
ergy quasiparticles have arbitrarily large lifetimes r k . On 
the other hand, these very long-lived quasiparticles, hav- 
ing low energies, contribute little to the energy current. 
Hence k/T is finite. 



INTERACTING CASE 

For T > 0, the dangerously irrelevant interaction V 
must be taken into account. The most important effect 
of interactions is to shift the phase transition, such that 
the QCP is approached at finite temperatures from the 
diffusive metal phase, as shown in Fig. [TJ This is cap- 
tured by a renormalized mass i?, discussed below, which 
is positive above the transition. R acts as an effective gap 
for low energy quasiparticles, thus rendering all dc trans- 
port coefficients finite. The situation is in contrast with 
Ref. Q , where interactions regularize transport by intro- 
ducing quasiparticle scattering. Here, the leading low T 
dc conductivities are obtained from a Hartree-Fock (HF) 
analysis, where the only effect of interactions is to shift 
the quasiparticle mass R. The HF results are shown in 
the last column of Table Q] These results are valid for 
extremely low temperatures, such that lnlnA/T ^> 1. 
To study transport on a much broader region surround- 
ing the QCP, we introduce a classical treatment of the 
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order parameter which, when supplemented by Langevin 
dynamics, will be shown to capture the correct quan- 
tum critical transport behavior over the region where the 
much weaker condition, In K/T ^> 1, is satisfied. 



Classical action for order parameter 



values of s — s c . This is due to the absence of long- 
range order (LRO) in 2d at finite T. For an 0(2) order 
parameter, however, quasi LRO is established at a T > 
Kosterlitz-Thouless transition. 

A description in terms of a classical action [12j is ap- 
propriate whenever log A^/T ^> 1. In this limit, U <S^T, 
so that modes with uj n ^ are significantly gapped. 



With interactions, the critical point is shifted away 
from s c = 0. To linear order in V, 
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In the vicinity of the QCP, we renormalize V by a one 
loop RG equation, up to the scale where the system either 
develops a gap, or when the rescaled temperature reaches 
an upper frequency cutoff A w [7. 11|. 
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The static properties of the finite T model can be stud- 
ied by integrating out all non-zero Matsubara frequency 
modes. After rescaling, <f> = \ / 2mTijj, r = x/V2m, the 
Lj n = mode has the following classical action, 



S r = I d 2 X 



where U = 2mTVn. This theory is super-renormalizable, 
and is rendered UV finite by introducing a renormalized 
mass R. 



R = R-2U 
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R has a universal expression in terms of s — s c , reflecting 
the contribution of the u>„ ^ modes 
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(15) 



where F(s — s c , rf) — > In y/l + rj 2 as s — > s c . Solving this 
self-consistent equation at s = s c yields, 



R 



(16) 



up to a prefactor of order lnln(A w /T). Note that, for 
T > 0, R is always positive, even for arbitrarily negative 



Dynamics of order parameter 

We approximate the low frequency dynamics of the 
classical order parameter by a Langevin equation (model 
A dynamics of Ref . [13( ) , 



dt 



(17) 



(f;(x,t)f n (x',t')) = 2 v S(x-x')S(t-t>) 



Equal time correlators computed with these dynamics 
are equal to those of the classical action S c , as necessary. 
The appearance of the "bare" value of rj in cq. (JTTJ) is 
due to the fact that dispersion in the quantum action is 
non-local in time and therefore is not renormalized. 

Consider the HF approximation, in which (fT"7|) be- 
comes a linear equation with mass R. Solving for cr, 



Recr(ft) = 




For R = 0, this reproduces the non-interacting result, as 
expected. On the other hand, for R ^ 0, we obtain a 
finite dc conductivity, 



m — 

T 



(18) 



We note that Eq. fTS]) disagrees with Ref. [14(, which 
predicts a ~ lnln(A w /T). This is due to an erroneous 
computation of T* in Eq. (12) of that reference. 

Naive use of Eq. (jTTJ) yields divergent values of a and 
k/T. This is not surprising: the Langevin equation as- 
sumes classical modes, whose occupation factors satisfy 
equipartition, / cq (w) = T/uj. However, inspection of the 
Boltzmann approach, eq. (|13p . shows that for such distri- 
bution, a and k/T have UV catastrophes. In this sense, 
the Langevin equation does not capture the correct dy- 
namics of high energy quantum modes. However, high 
energy modes are very weakly perturbed by the quar- 
tic interaction. Thus, the correct result is given by the 
Boltzmann equation with a full Bose distribution, as in 
Eq. (fT3|) , but with a chemical potential set by R. Equiva- 
lently, this corresponds to evaluating the one-loop quan- 
tum expression ([5]) with chemical potential R. This yields 
the last column of Table Q] 
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The HF results can be obtained independently from an 
exact solution of the quantum model ([2]) in the large N 
limit [l5| , where N is the number of components of the 
order parameter (N = 2 for superconductivity). This 
is an important check that the Langevin equation |T7|) 
captures the correct universal dynamics. We see from 
Eq. (fl"5|) that, for s — s c and at low T, 



U 1 
R * WrjR 



1 



lnlnA w /T 



which justifies HF provided that lnln(A u /T) is a large 
number. More generally, to go beyond HF, we must con- 
sider higher order corrections in U . From the Kubo for- 
mula and the fluctuation-dissipation theorem, we deduce 
that the dc electric conductivity obeys 



(2e) 2 T 
h 8ttt]R 



(19) 



for some scaling function Q a satisfying $ CT (Q,ry) = 1. A 
similar analysis applies to a and fc/T, with the important 
difference that substractions are necessary to cancel lead- 
ing UV divergences, as discussed above. When working 
with the renormalized R, the only UV divergence comes 
from the diagrams already computed. Thus, 



2e (U 

(1 — ^quantum, 1 loop T" ^ a 1 J V 
K> ^quantum, 1 loop . <J) ( 



where the small R limit of the (quantum, 1 loop) results 
§^jln| for a, and Eq. © for k/T. The functions 



is 




$ 7 can be evaluated numerically by introducing a lattice, 



U/R 



FIG. 2: Estimated scaling function <& a {U/R,ri = 1) from 
numerical integration of the Langevin equation (|17[) . Con- 
vergence for static quantities of the Langevin algorithm was 
tested by comparison with results from a Wolff cluster algo- 
rithm. Each data point represents at least 2 x 10 7 time steps 
on a lattice of spacing a — 0.1 and linear size L — 64, which 
for static quantities with n = is known to approximate well 
the continuum limit a — > 0. Results were normalized by the 
U/R = result and fit by a quadratic polynomial (solid line) . 
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and requiring that the renormalized R be the same in the 
lattice and continuum theories, 



Ri 



R-2U 



2tt 



2-7T 4 — 2 cos k T — 2 cos k,. 
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Fig. [5] shows the scaling function $> a (U/R,r] = 1). 
This, combined with Eqs. (fTB")) and p^|) . gives the elec- 
tric conductivity for the entire quantum critical regime. 
We stress that these results rely only on the condition 
logA/T 3> 1. This use of Langevin dynamics to obtain 
full scaling functions for transport quantities at a QCP 
should be applicable at many other transitions. Here, we 
have used them to find WF violation and an anomalous 
thermoelectric conductivity at a transition of experimen- 
tal interest. 
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